Unitarity in periodic potentials: a renormalization group analysis 
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We explore the universal properties of interacting fermionic lattice systems, mostly focusing on 
the development of pairing correlations from attractive interactions. Using renormalization group 
we identify a large number of fixed points and show that they correspond to resonant scattering in 
multiple channels. Pairing resonances in finite-density band insulators occur between quasiparticles 
and quasiholes living at different symmetry-related wavevectors in the Brillouin zone. This allows a 
BCS-BEC crossover interpretation of both Cooper and particle-hole pairing. We show that in two 
dimensions the run-away flows of relevant attractive interactions lead to charged-boson-dominated 
low energy dynamics in the insulating states, and superfluid transitions in bosonic mean-field or XY 
universality classes. Analogous phenomena in higher dimensions are restricted to the strong coupling 
limit, while at weak couplings the transition is in the pair-breaking BCS class. The models discussed 
here can be realized with ultra-cold gases of alkali atoms tuned to a broad Feshbach resonance in an 
optical lattice, enabling experimental studies of pairing correlations in insulators, especially in their 
universal regimes. In turn, these simple and tractable models capture the emergence of fluctuation- 
driven superconducting transitions in fermionic systems, which is of interest in the context of high 
temperature superconductors. 



I. INTRODUCTION 

Fermionic ultra-cold atoms with nearly resonant scat- 
tering in the unitarity regimei"— realize the strongest 
possible form of Cooper pairing, revealed by criti- 
cal velocity^^ii^. It is natural to expect that zero- 
temperature normal states near unitarity would be very 
strongly correlated with a host of unconventional proper- 
ties, possibly bearing some resemblance to those found in 
cupratesi^. The unitarity limit is therefore an excellent 
starting point for studies of correlated fermionic super- 
conductors and insulators, which has not been exploited 
enough in literature. The benefits are both theoretical 
and experimental. Systematic perturbative and renor- 
malization group calculations are feasible mainly because 
the unperturbed ground state (fixed point) is a simple 
state, the vacuum or a band insulator. Experimentally, 
the unitarity limit is routinely accessed in cold gases of 
alkali atoms tuned near a broad Feshbach resonancei^ii^. 

Our ultimate goal is to address the long-standing ques- 
tions about the nature of unconventional normal states 
proximate to strongly paired fermionic superfluids or su- 
perconductors. We design here a simple and tractable 
model in which the fermionic excitation gap is opened by 
an external periodic potential, rather than strong inter- 
actions. Attractive interactions between quasiparticles 
whose energy scale exceeds this gap can still give rise 
to pairing and superfluidity. The need for strong inter- 
actions justifies asking if the normal state proximate to 
the superfluid might have some unconventional proper- 
ties reflecting strong correlations, especially in universal 
regimes shaped by resonant scattering. 

The inquiry into unconventional superfluidity from the 
resonant scattering point of view began a long time agc^^. 
The present interest in this subject is driven in parallel by 
a variety of unconventional superconductors in condensed 



matter physics, and ultra-cold atoms with nearly reso- 
nant scattering. The recent theoretical studies of scatter- 
ing resonances in lattice potentials'^"— often rely on two- 
channel tight-binding models, featuring fermionic atoms 
resonantly coupled to closed channel bosonic particles. It 
has been argued that models of this kind provide a good 
effective description of the microscopic lattice systems of 
interes1i2^"— . The flndings of these studies include lat- 
tice Feshbach resonances shifted from their empty-space 
values. 

Most of the mentioned theoretical works approach res- 
onant scattering from a somewhat microscopic angle, ex- 
emplifled by perturbation theory with the vacuum un- 
perturbed ground-state. Indeed, the usual universal be- 
havior of particles tuned to a broad Feshbach resonance 
is established in the low density limit. However, uni- 
versality is a many-body phenomenon and its complete 
description requires field-theoretical tools such as renor- 
malization group. This issue becomes pressing in the 
present problem of interest, a band insulator whose non- 
interacting ground-state contains a macroscopic number 
of fermions in fully populated bands. These fermions 
may not be dynamically inert despite the Pauli exclusion 
principle, due to the strong interactions which bring the 
system to its unitarity regime. 

In this paper we take a field-theoretical approach to 
nearly resonant pairing between gapped fermions. Aban- 
doning all microscopic details in a renormalization group 
(RG) calculation allows us to gain a perspective on the 
generic and universal behavior of a large class of fermionic 
lattice systems. The effective theories we subject to 
RG are constructed to preserve the universality class of 
the microscopic system. This ensures that the univer- 
sal phase diagram and other macroscopic properties of 
the microscopic system are correctly captured despite 
the neglect of microscopic details. The price to pay is 
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the necessity to deal with multiple flavors of low energy 
quasiparticles, such as particle and hole excitations which 
may exist at multiple wavevectors in the first Brillouin 
zone. Note that the two-channel models with dynami- 
cal bosonic fields used in many previous studies are not 
guarantied to be in the same universality class as the mi- 
croscopic system, and may in some cases describe differ- 
ent physics than this paper. We characterize the univer- 
sality stemming from resonant scattering of quasiparti- 
cles in band insulators, and discover generalized unitarity 
regimes in which quasiparticles of different flavors scat- 
ter resonantly. The manifestations of unitarity which we 
discuss include universal ratios of measurable quantities 
such as critical temperature, pressure and density. We 
also analyze the types and conditions for pairing instabil- 
ities, conventional versus unconventional superfluid tran- 
sitions, and emphasize the existence of correlated bosonic 
Mott insulating states in the phase diagram. 

The RG analysis reveals why pairing fluctuations in- 
deed play the crucial role in systems of gapped fermions 
with short-range attractive interactions. Unlike previous 
studies of unitarity in continuum lattice potentials, which 
focused on the zero-density limit^Si"— , we point out that 
the unitarity regime in the same universality class can 
be found at flnite densities near any zero temperature 
band insulator to superfluid transition. The structure of 
fixed points depends on whether both particles and holes 
participate equally in the dynamics, or just one of the 
two quasiparticle types. In the latter case, the exact RG 
equations can be derived, which allows one to track the 
run-away fiows of attractive interaction couplings (as in 
some studies of Iron-pnictides^i^) . It is these run-away 
flows that can lead to boson-dominated dynamics at low 
energies. We find that instabilities in the particle-hole 
channel are discouraged by attractive interactions. 

The run-away flows imply the ultimate RG breakdown 
when the diverging couplings reach cut-off scales. How- 
ever, the resulting low energy bosonic dynamics is known 
to introduce additional flxed points associated with su- 
perfluid transitionsSi, which appear as strong-coupling 
flxed points in the present RG. The superfluid transition 
in this regime can be either in the bosonic mean-field, or 
XY universality class. The mean-field universality with 
dynamical exponent z = 2 emerges as the result of run- 
away fiows from the unitarity dominated by either parti- 
cles or holes, while the z = 1 XY universality is related 
in the same fashion to the unitarity shaped by both par- 
ticles and holes. Therefore, the analysis here provides a 
glimpse of the more complete structure of fixed points 
in theories of fermionic particles with attractive interac- 
tions, sketched in FigHJ 

The mentioned finite-density fixed points describe uni- 
tarity in zero-density effective theories of particle and 
hole excitations. Therefore, one can relate a nearly crit- 
ical interaction strength to scattering lengths in colli- 
sions among particles and holes. Any attractive inter- 
action in two dimensions effectively puts the low energy 
quasiparticles into their Bose-Einstein condensate (BEG) 



u 




FIG. 1: (color online) A hypothetical renormalization group 
(RG) flow diagram at zero temperature of a fermionic lat- 
tice theory with attractive interactions in d > 2 dimensions. 
The parameters are fermion density-density interaction U and 
bandgap (or negative chemical potential) Eg. This paper ex- 
plores in detail the vicinity of two weak-coupling fixed points 
which govern the pair-breaking superfluid to insulator tran- 
sition: Gaussian (G) and unitarity (U), separated in propor- 
tion to e = d — 2. The negative couplings U below unitarity 
experience run-away flows under RG and lead to the boson- 
dominated dynamics. A superfluid-insulator transition in this 
regime is captured by a bosonic eiTective theory with dynam- 
ical exponent z = 2 (mean-fleld universality) or z = 1 (XY 
universality). In the latter case, one expects an additional 
fixed point (XY) in d < 3. The Gaussian fixed point of the 
bosonic eiTective theory appears at the transition line in the 
limit U — > — oo, Eg — >■ oo. The shaded area is the superfluid 
or superconducting phase, and the red thick line is the second 
order superfluid-insulator transition. The dashed green line 
encloses the region in which the fermionic RG is valid. 



limit, so quasiparticles injected in the insulating state 
immediately combine into bound-state pairs^^ (whose 
size can be very large at weak couplings). The effec- 
tive Bardeen-Cooper-Schrieffer (BCS) regime exists only 
above two dimensions, at least in the weak coupling limit 
{\U\ < \U* \ (xe in Figd]). Note that localization tenden- 
cies due to the lattice potential enhance the strength of 
effective interactions with respect to those between com- 
pletely free fermionsl^'-''. 

We begin the discussion by laying out the effective the- 
ory of a lattice fermionic system in the section IF^ Then, 
in the section Hi Al we analyze the exact RG for a single 
species of fermions, either particles or holes, and reveal 
the development of short range pairing correlations in in- 
sulating states. The section lTl Bl m'csents the RG calcula- 
tions for particles and holes and identifies a large number 
of flxed points associated with resonant scattering. This 
analysis is expanded in the section lll CI to multiple quasi- 
particle species living at different wavevectors in the Bril- 
louin zone. All results and conclusions are summarized 
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in the discussion section Hill 



A. Model 

As a generic model of a band insulator we consider the 
imaginary-time action of neutral fermionic particles with 
interactions C/, in a lattice potential V{r): 
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where the summation over spins a S {t;i} is implicit. 
This is a microscopic multi-band model defined in contin- 
uum space and not a priori tied to the vicinity of a critical 
point. The density of particles is tuned to any number of 
completely populated bands at zero temperature by plac- 
ing the chemical potential in a bandgap. The dynamics 
of the resulting band insulator can be described by an ef- 
fective theory of low-energy quasiparticles belonging to 
the valence and conduction bands. Sufficiently strong 
attractive interactions U can drive the system into a su- 
perfluid state, and a similar instability can be created 
by bringing the chemical potential sufficiently close to 
a band edge, even at weak couplings. A qualitative ex- 
ample of the superfluid-insulator transitions is shown in 
FigH 

In order to derive the effective theory we formally in- 
tegrate out the fermion fields from high-energy bands 
in the path integral. At best, this can be done pertur- 
batively, for example using the Feynman diagram tech- 
nique. Unless the perturbative integration breaks down, 
the effective theory takes form 
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in terms of the quasiparticle fermion fields /. In this 
paper we will consider cases in which the band indices 
riijiTii denote one or two bands. The most generic tran- 
sitions are driven by the chemical potential fj,, so that 
only one band is important (see Figl^l (p) and (h) tra- 
jectories). The transitions involving particles and holes 
(the (ph) trajectory in FigH]) require at least two bands 
for a complete description. 

A potential problem recognized in a number of stud- 
ies is that in the vicinity of resonant scattering the mi- 
croscopic interactions may correspond to energy scales 
(much) larger than the bandgap. Then, many high en- 
ergy bands may be significantly hybridized with the con- 
duction and valence bands. In the present formulation of 
the problem, this can lead to a strong renormalization of 
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FIG. 2: (color online) Superfluid transitions (thick red lines) 
out of a two-dimensional band insulator at zero temperature. 
The three shown transitions correspond to arbitrarily chosen 
different strengths of contact attractive interactions, becom- 
ing stronger going from top to bottom. The lattice potential is 
given by V{r) — 2V[cos{2-kx /aL) -I- cos(27ry/aL)], where V is 
the lattice amplitude and ul the lattice spacing. Both V and 
^ are measured in the units of "recoil energy" 
The thin red line outlines the band edge of the corresponding 
non-interacting model. The dashed blue lines are trajectories 
in the parameter space along which the transitions dominated 
by particles (p), holes (h), or both (ph) can occur. 



low-energy quasiparticle dispersions En{k) and effective 
interactions U™]^'^ . In most circumstances this does not 
endanger the analysis which follows. 

The goal of this paper is to explore the universal prop- 
erties of nearly resonantly interacting fermions in peri- 
odic potentials, for the purpose of which we apply renor- 
malization group. Therefore, the precise functional forms 
and values of the effective En{k) and U™]^"^ are not 
of interest, and we do not attempt to derive them from 
any microscopic model. The perturbative expansions of 
En{k) and in powers of the ratio U/Ehe between 

the microscopic interaction strength U and the small- 
est energy Ehe of the integrated high-energy fermions, 
need not converge fast. However, these expansions must 
be convergent, otherwise the effective En{k) and U™]^^ 
would contain singular features which would invalidate 
our analysis. In general there should be a minimum num- 
ber of bands which must be kept in the effective theory 
in order for it to be properly analytic, and this number 
may grow when one approaches the resonant scattering 
in empty space at V{r) = 0. We shall assume that this 
number is never larger than two in the present cases of 
interest. An assumption of this kind is made in all other 
studies and it is justified by the fact that a lattice poten- 
tial shifts the scattering resonance from its empty-space 
position toward the effective lattice BEC limitii^i^i. In 
other words, the effective unitarity limit in a lattice which 
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we seek to describe corresponds to the BCS regime in 
empty space for the same microscopic interactions, where 
U/Ehe is not too large. 

II. RENORMALIZATION GROUP ANALYSIS 

Here we apply renormalization group (RG) to a band 
insulator in the unitarity regime. We will identify various 
fixed points associated with unitarity (resonant two-body 
scattering) which emerge in the presence of an external 
periodic potential V{r), but otherwise are analogous to 
the unitarity fixed point of a uniform system at V{r) = 0. 
The main difference is that the fixed points we shall dis- 
cuss occur at finite densities of microscopic particles, cor- 
responding to fully occupied bands, whereas universality 
in the uniform system stems from a zero density fixed 
point. 

There are two characteristic situations which will be 
considered separately. First, one fermion species (ei- 
ther particles or holes) generally dominates dynamics in 
the unitarity regime, so the renormalization group equa- 
tions can be derived exactly to all orders of perturba- 
tion theory. This is extremely useful because run-away 
flows of interaction couplings can be traced more reliably. 
The second situation is more special and occurs when 
both particles and holes participate equally in dynam- 
ics. Then, the fixed point structure becomes intricate, 
but can be accessed only in an e (or large- TV) expansion. 
At the end we briefly discuss extensions to more realistic 
cases with multiple relevant fermion species, and mea- 
surable manifestations of the universality class. 

A. Transitions involving one fermion species 

A transition dominated by either particles or holes, 
but not both, is generally caused by chemical potential 
changes as illustrated in Figl2] with (p) and (h) dashed 
lines. As a natural starting point one can imagine a band 
insulator either in the deep BCS limit, or with a very deep 
lattice potential, where the chemical potential is brought 
much closer to one of the bands than to the other. How- 
ever, such extreme regimes are not necessary initially be- 
cause they are created by the RG flow (the bandgap is 
a relevant operator). We shall discover that an attrac- 
tive interaction undergoes a run-away flow in two dimen- 
sions and competes with the bandgap. This competition 
is resolved at cut-off scales where the RG breaks down. 
However, a new strongly-coupled universality class takes 
over in that limit, associated with superfluid to Mott- 
insulator transition with dynamical exponent 2: = 2 at 
the intersection of the thick red and dashed-blue lines 
(p) or (h) in FigEl 

A characteristic weak-coupling interacting fixed point 
contained in a theory of interacting fermions is unitarity, 
and it is found at zero temperature when the chemical 
potential lies exactly at the boundary between a fermion 
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FIG. 3: The RG flow of the interaction coupling U in the 
theory Q. 

band and a bandgap (or vacuum). The effective action 
contains a single species of interacting spinful fermions, 
not much different from the theory of the system with- 
out a lattice. If these fermions live in a valence band, 
we can immediately reformulate the theory in terms of 
holes. Therefore, we can always write the critical effective 
theory in the limit of zero quasiparticle density. The uni- 
versality class at unitarity will be the same as discussed 
in Ref.36.37. 

The critical theory of interest in d dimensions, aug- 
mented by the relevant chemical potential is: 

^1 = / mfl^i-^u: + Eik))f,^^ (3) 

where k = (w, fc), Vk = da;d''fc/(27r)''+\ and 

E{k) = Eo + -^. (4) 
Zm 

All fermion loop Feynman diagrams vanish (the poles of 
the Green's functions on the loop lie in the same complex 
half- plane), and only the interaction coupling is renor- 
malized by a summable geometric progression of ladder 
diagrams. The exact RG equations are found to be^: 

^ = 2E, , ^ = (2-d)C/-^C/^ (5) 

where I is the scale parameter. Eg = Eq + U is the effec- 
tive bandgap, and H is a positive cutoff-dependent con- 
stant. A fixed point is always found at Eg = 0, U — 0. 
An additional non-trivial fixed point is found in d ^ 2 at 
Eg ~ 0, U — U* = (2 — d)H^^, which describes attractive 
interactions if d > 2. The schematic flow of interaction 
couplings is shown in FigjHl 

Any attractive interaction U < in d ^ 2 has a run- 
away flow to U — )■ —00, while in d > 2 it needs to be 
large enough to flow toward U — >■ —00. Repulsive inter- 
actions U > 0, on the other hand, flow to the Gaussian 
fixed point in c? > 2. Therefore, only attractive interac- 
tions can produce strongly correlated states. Since the 
RG equations are exact, we can precisely characterize the 
run-away flow, assuming that the effect of multi-body col- 
lisions remains negligible at least in the insulating state. 
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Solving for U{1) in two and three dimensions we obtain: 

^ ' ] u(o) d - 3 I 

In both cases, the run-away flows have vertical asymp- 
totes so that [/(/) diverges at a finite value oi I {I — 
|n[/(0)|~^ in c? = 2). This indicates that Cooper pairs 
become stable at a finite length scale and have a finite 
coherence length despite the fermion bandgap. However, 
the interpretation of the run-away fiow breaks-down at 
the cut-off scale because so large U will pair-up the high- 
energy fermions, which were assumed to be unpaired in 
this RG procedure. A boson-dominated dynamics takes 
over the shortest length-scales under consideration. 

Note that U{1) — >■ — oo at a finite I cannot be immedi- 
ately interpreted as a signal of superfiuidity. This is due 
to the fact that at finite I we do not yet have a theory 
which transparently describes dynamics at macroscopic 
scales, while superfiuidity is verified only by macroscopic 
long-range correlations. Since RG is based on integrat- 
ing out high energy modes, it does not provide a precise 
answer to the question of what phase the system lives in, 
but only gives an indication. 

The fermion gap in an insulating state grows exponen- 
tially under RG, Eg{l) = Eg{0)e^'. If Eg{l) is the first to 
reach the cutoff scale, the further RG flow is halted (RG 
breaks down) in a state apparently devoid of particles. 
This is a band insulator. It is obtained in 2+e dimensions 
for any sufficiently weak interaction |C/(0)| < \U*\ ~ e, 
or even for \U{0)\ > \U*\ provided that the gap Eg{0) 
is large enough. If U{1) is the first to reach its cut- 
off instead, then boson-dominated dynamics at shortest 
length-scales requires switching to a purely bosonic effec- 
tive theory in order to determine what happens at large 
length-scales. Both insulating and superfluid phases are 
possible in this limit, despite a finite fermion gap, but 
the transition between them is in a different universality 
class than the BCS pair-breaking transition. 

B. Transitions involving particles and holes 

A special case is obtained in the vicinity of vanish- 
ing gaps for both particle and hole excitations. A pair- 
ing transition influenced by a corresponding flxed point 
cannot be obtained by changing the chemical potential 
alone, but can be accessed by tuning interaction strength 
or lattice depth, at a fixed particle density (see FiglJl 
the (ph) dashed line). The values of /i and V should 
lie at the intersection of effective conduction and valence 
bands where the effective bandgap closes. We shall again 
discover a run-away flow of interaction couplings, but 
this time it quickly invalidates the perturbative RG. The 
run-away flow is expected to eventually lead to a strong- 
coupling flxed point in the XY universality class, asso- 
ciated with the superconductor to Mott-insulator transi- 
tion at an integer number of bosons per lattice site. 



The critical theory for valence (v) and conduction (c) 
electrons is: 

S2 = Y.h^ fUo. + En{k)) (7) 

X f'f f f 

7ni, ki-\-q,a'' ni,ki,a'' 1712,^2— q,i3-' 712, k2,^ ' 

where n € {c, v} and 

E4k)^-E,Q-- , E,{k) = E,o + ■ 

2mv 2mc 

Here we made the simplest assumption that the bandgap 
Eq = Evo + Eco > is direct and small (or vanishing) at 
only one wavevector in the Brillouin zone. Modifications 
of this assumption are straight-forward and the appro- 
priate more realistic circumstances will be discussed in 
the following section. Performing a particle-hole trans- 
formation for the valence band cannot help us construct 
an exact RG procedure. Instead, it is convenient to work 
directly with the native particle degrees of freedom. 

The interaction couplings C/^J^™^ in the band repre- 
sentation are derived from microscopic short-range in- 
teractions in real-space. For example, a pure contact 
interaction Ui(jl^{r)t(;^(r)ipj^{r)ip i^{r) in ([T]) gives: 

c/™^7(fci,fc2,g)= (8) 

U J d'^r w;;,,,fci+q(r-)w„^_fc^(rX,^ fc^„g(r)u„^ fe^(r) 

uc 

where UC indicates integration over the lattice unit-cell, 
and il)n,k{i') = Un,k{i')e.^^^ are Bloch wave-functions. 
This expression illustrates an important property of in- 
teractions in the band representation which follows from 
the overlap features of the Bloch wavefunctions. As a rule 
of thumb, the couplings J7^^™^ are largest by magnitude 
if tii = rrii for both i = 1, 2 and smallest if Ui ^ for 
both i — 1,2. The strongest interaction channels involve 
a single band, while the interband couplings are weaker. 
This is a natural situation for generic band structures 
and short-range interactions, but it could be reversed in 
principle. 

It is important to note that the spatial dependence 
of the microscopic interaction potential U (r) on the dis- 
tance r between the interacting particles is not automat- 
ically irrelevant (in the RG sense) in the presence of the 
lattice. Short-range variations of J7(r), at or below the 
lattice spacing length-scale, affect the relative strength 
of the couplings U™]^^, which may lead to non-trivial 
interacting fixed points as discussed below. Only the 
long-range variations and the related crystal momentum 
dependence of the interaction couplings are not relevant 
in the vicinity of the flxed points of interest. 

After normal ordering the interactions take the form: 
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The generated quadratic terms /7™ = no/2 x [/^^ , 
where no is the particle density in the ground-state, 
are similar to a "charging energy" since they effec- 
tively shift the chemical potential as a result of inter- 
actions. However, they also couple the particles in the 
two bands, so we must redefine the bare Green's func- 
tion. One way to accomplish this is to treat as a self- 
energy correction to the non-interacting Green's function 

[iuj - En{k)]^'^Snm- 



iuj-Ec{k)-U^ -Ul 
[iLo — Zi{k)){iuj — Z2{k)) 



(10) 



It is convenient to define C,n{k) — En{k) + and ^ 
y^U^ (note that = {U^)*). Then: 



2i/2(fc) 



Cc + Cv 



± 



Cc-Cv 



(11) 



are the new poles of the bare fermion excitations, and 

5;r(fc, lUj) = (iLO - Cc - Cv + Cn)Snm + ^(1 - <5„™) • (12) 

It is easy to show that both poles are always real, one 
being positive (particle- like) and the other negative (hole- 
like) as long as > CcCv- We will assume that the 
latter condition is satisfied, so that the system is a band 
insulator despite the "charging energy". Consequently, 
we can expand the poles up to 0{k'^): 
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where are the bare quasiparticle gaps, and Mi are the 
quasiparticle masses given by: 



M^^2 = a X 
The parameter 
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captures the amount of mixing between the two bands 
(0 < a < 1); £'g = (Cc — Cv)|fc_o effective fermion 

bandgap. For a — 1 there is no band mixing since Mi g 



{TOcTTiv}. In general, a > |/3| = \mc — mv\/{mc + m^) is 
required in order for both Mi to remain positive. Other- 
wise, band inversion is caused by large interband cou- 
plings and it must be taken into account by redefining 
the low-energy quasiparticles, which then live at some 
different momenta in the first Brillouin zone. We shall 
come back to this situation at the end. 

Now we set up the RG. As usual, we keep the masses 
TOc and TOv in the absence of interactions fixed under RG. 
While this does not imply that Mi will be fixed, it sets 
the scaling dimension for the field operators to d/2. The 
scaling of coordinates and couplings 



r' = (16) 



is followed by the diagrammatic integration of high- 
energy fields living at all Matsubara frequencies and mo- 
menta within a shell |fe| G (Ae^'^', A), where A is a cut-off 
momentum scale and dl is an infinitesimal increment of 
the scale parameter /. The resulting one- loop renormal- 
ization of the quadratic and quartic couplings is summa- 
rized in table m The relevant cutoff-dependent renormal- 
ization scales are: 
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The RG equations involving all four f/™ and all sixteen 
^nin™^ couplings (uot all of which are independent) are 
too complicated to be fully solved. A part of the problem 
is that the parameter a can also flow under RG, as a 
result of the renormalization of the couplings ?7™. In 
order to simplify notation let us absorb the bare fermion 
gaps into the quadratic "charging" couplings since they 
flow the same way under RG: C/" — ^ C/" -I- i5„o- We 
begin by noting that the RG equation for the quadratic 
couplings is: 



= 2[/; 
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e dimensions the interacting fixed points will 
oc e, implying [7™ cx e. Finite values for 
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T^'-^feni '^n2i '-^feni '-^ ln2 



Ery-/ (rjmiljjkm2 , jjinii-iljjm2k 
^2kl[^LJ„j^k^n2t +^nik'^ln2 



,rrlm-i rrkm2 
''-^kni '-^n2l 



fcni Zn2 



j_Tjlm.2rrkmi , rjlm2Tr 



mik 
Ini 




c\ Mm >v cX Mee'" *^ X Mee 





TABLE II: Relevant interaction vertices near the zero- 
quasiparticle-density fixed points for a = 1 (defined in the 
text). 



plings U^^ and U^^ are larger than the interband U^^ 
and U^^, so the attractive intraband channels dominate 
at macroscopic scales and lead to Cooper pairing even if 
the interband channels are repulsive. Instabilities in the 
particle-hole channel are possible only if all interactions 
are repulsive, or if for some reason U^^ is repulsive and 
stronger than the attractive intraband interactions. 

In d = 2 + e dimensions with e > it is convenient to 
define: 



Ery-r (rjmiljjkm2 , j- rmi / r T-m2 fc 
.rrlmirjkm2 I rjlm-iTjm2k 



m(l + /3) , m^ = m{l-l3) , /3 = 



and the rescaled independent dimensionless couplings 

(Wc; '^V; '^CV5 ; '^C; ^g)- 



El^n ( J jmxra2 T jkl . TTmim2Trlk 
kl ^ 

m2^1Tjlk \ 



.Trm2mxrrKl \rr 
'^'-'ik '^nin2 T" I 



TABLE I: One-loop diagrams which renormalize the couplings 
Un (first two) and U^^\^^ (last four). The renormalization 



constants K-\k, K'oi 



and K^kk' 



are given m 



all f/™ in d > 2 dimensions uniquely determine the value 
for a, which has to be fed back into the RG equations 
to self-consistently determine the fixed points. This can 
be done only numerically. However, analytical solutions 
for a subset of fixed points can be found if the couplings 
U^, t^cc J f^cvj f^vv: ^cv 8.nd their complex conjugates are 
all zero. In this case, it follows from (jisp and ([T8| that 
a — 1 and does not flow under RG. This will be the focus 
of the following discussion. 

In two dimensions there is only one weak-coupling fixed 
point, at U^' = 0, C/™i„™^ = 0. The flow of U^^, C/;; 
and U^^ — U^^ is of the same type as shown in FiglH 
at d = 2, so that the attractive interactions undergo a 
run-away flow, while repulsive interactions flow to zero. 
The interband interaction U^^ = U^^ has the opposite 
behavior, a repulsive one keeps growing, while an attrac- 
tive one flows to zero. In normal circumstances, due to 
the Bloch wavefunction properties, the intraband cou- 



= Ke 



^ C V 
Ur 



1 + 



(19) 



Ke 



1-/32 

i9 



JJVC _^ JJC 



Ke- 



1-/32 



-ie 



U^^, = Ke 



where K = [2itY / [Sd^'^rn). These interactions are dia- 
grammatically represented in the Table HIl The RG equa- 
tions for a = 1 are: 



dwc 
"dT 
duv 
"dT 

dUcv 

"dT 
d; 

dwc 

"dT 

deg 
dl 



= 2e„ 



-Uc - Aul - Aul 



-Mcv + 2ul, + 8(1 - P'')ul 

-M,„ - + 4Ucv1ti] 

2liv 2l/rv 



(20) 



4Mn 



1-/3 1-/32 I_p2 



Above two dimensions, there are seventeen flxed points 
with a = 1. Sixteen of these flxed points Fi — Fiq are 
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given by all possible combinations of: 



tMcv>0 



tMm<0 



Uc e 



„„,,,{(>,o).(l.i),(o,4),(o,o)} 



Note that here is always zero. The RG eigenval- 
ues in the subspace of {uc,Uv,Ucv,u^) are ±e at all 
of these fixed points, and allow enumerating Fi — Fiq 
simply by the variations of relevant /irrelevant flows of 
(uc, Mv, Wcv, Wm)- This is illustrated in Fig |4] for the 
first fifteen fixed points at which at least one of the 
Uc, Uv, Ucv, couplings is zero. Only the Gaus- 
sian fixed point is fully stable, while Fiq with all Uc, 
Mv, Ucv, non-zero is fully unstable. The cou- 
pling Uc is irrelevant only at the Gaussian fixed point, 
while it is found to be marginal at (uc, Mv, Wcv, "m) e 
{(-1/4, 0,0,0), (0,-1/4, 0,0), (0,0, 0,-1/4)} and rele- 
vant otherwise. 

The remaining fixed point Fn is the only one with 
Uc > 0: 
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64 ( 11 - 4/32 + 8^4/3'' -7/32-^4) 
1 3 



, Um = 



Uq — — ^cv 



8(5- 413^ + 2^/4^4 - 7/32 + 4^ 

It has only one relevant direction with RG eigenvalue e, 
the Ue component being the largest in the corresponding 
eigenvector. The RG flow in the vicinity of this fixed 
point is illustrated in Fig El Note that formally there are 
other solutions stemming from ([20)1 . but they have u^ < 
corresponding to time-reversal symmetry violations. 

Whenever Uc, u^, Ucvi relevant, their RG eigen- 

value is e, the same as the RG eigenvalue at the unitarity 
fixed point of a uniform zero-density system which corre- 
sponds to vacuum resonant scattering. While this is not 
surprising when it comes to pairing of two quasiparticles 
Uc or two quasiholes Uv, it is interesting to note that the 
same resonant scattering interpretation can be applied to 
the couplings Ucv and Mm. We can identify the resonantly 
scattering quasiparticles by tuning to a fixed point with 
only one finite coupling and then taking a closer look at 
the operator corresponding to that coupling. In the case 
of Ufii, the operator is 



J caJ vl3 J cj3 •> va •/ caJ vl3 ■> vaJ c(3 



to I 



where the operator $s ^ (/cf/vj, - /c;/vt)/%/2 annihi- 
lates an interband singlet and the operators = /cf/vti 
^ti = fcifn and $to = (/ct/vi + /ci/vt)/\/2 annihilate 
triplet pairs. The fixed point(s) at Um < can now be 
associated with the resonant scattering in the interband 
singlet Cooper channel. Note that the absence of fixed 
points at itm > rules out resonant scattering in the 




FIG. 4: The fixed points and RG flow of interaction couplings 

(Mc, Uv, -iicv, Wm) at Me — 0: (a) Um = 0, (b) Ucv = 0, (c) 

Uv — 0, (d) Mc = 0. In front of the exposed "cube" faces 
are the run-away regions, signifying pairing correlations in 
various channels and possible symmetry-broken phases. If one 
extends the exposed "cube" faces in the directions of Uc > 0, 
iiv > 0, Ucv < and Um > 0, the obtained semi-infinite 
surface (whose corner is a shown "cube") encloses the basin 
of attraction of the Gaussian flxed point. 




Mcv>0 



FIG. 5: The fixed points and RG flow involving Ue 7^ for 
Um ~ 0, Uc = Uv The shaded semi-infinite surface encloses 
the basin of attraction of the Gaussian flxed point. 



attractive triplet channel (the fixed point at Ucv = 1/2, 
Um — 1/4 is fully repulsive in the particle-particle chan- 
nel). 

The interaction Ucv > at its resonant-scattering fixed 
point is repulsive in the particle-particle channel and can- 
not lead to a Cooper pair resonance. However, it becomes 
attractive in the particle-hole channel. Keeping only Ucv 
finite allows performing a particle-hole transformation in 
the valence band, after which the theory contains two 
similarly dispersing fermion fields (particles and holes) 
in their vacuum states, interacting attractively. Denoting 
the particle and hole annihilation operators as /„ = /^.^ 
and /j, = respectively, where a is the opposite spin 
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of a, the Ucv operator can be written as 

/ca/v^/v/j/ca ~ ~ faljjl pfa 

Now the operators Bs =_(/tA -_hh)/V^^ = /t/t' 
Bti = fifi and Bto = (/t/^ + /4,/f )/%/2 annihilate sin- 
glet and triplet particle-hole pairs. This interaction does 
not make any distinction between different spins, so that 
a scattering resonance appears simultaneously in the sin- 
glet and all triplet channels. The bound-state resulting 
from this resonance is an exciton, and symmetry breaking 
at finite particle and hole density can be either a singlet 
exciton condensate, or a ferromagnetic state depending 
on the other couplings as well as higher order terms in 
the action. In both cases, the present effective theory 
would favor ordering at zero wave vector, but the circum- 
stances discussed in the following section could lead to 
antiferromagnetic and other kinds of ordering at finite 
wave vectors. 

The behavior of does not fit this generic resonant 
scattering picture. Only at the Fn fixed point we find 
the flow of Me reminiscent of resonant scattering. The 
absence of other similar fixed points with 7^ 0, and 
the fact that the relevant direction at Fn is an almost 
even-amplitude linear combination of multiple couplings, 
indicate different physics: an "assisted scattering reso- 
nance" in the Cooper channel between a pair of fermions 
dynamically resonating between the conduction and va- 
lence bands. In fact, assuming — Qin a sufficiently 
strong interaction of this type would give rise to an ex- 
tended "sign-changing" s-wave superfiuidity in which the 
pairing gap on the conduction and valence bands has op- 
posite signs. Such an pairing is proposed to occur in 
iron pnictides^. Other kinds of pairings with different 
relative phases between the conduction and valence band 
pairing gaps could be obtained for other values of 9. 

The run-away flows in the vicinity of these fixed points 
are also very important. They indicate the kinds of in- 
stabilities of interacting fermions in lattice potentials and 
circumstances in which they can develop. This informa- 
tion has greater practical use than the detailed properties 
of the fixed points, because realistic systems can hardly 
be tuned very close to these fixed points (except cold 
atom systems which are tunable to the Uc = —1/4 and/or 
Mv = —1/4 fixed points). In generic lattice systems with 
attractive interactions we find that the favored phases are 
featureless insulators and superconductors. A singlet su- 
perconductor is indicated by the fiow of interaction cou- 
plings Uc, Uv and toward —00, although as emphasized 
in the previous section such run-away flows can also pro- 
duce bosonic Mott insulators in certain cases. Instabili- 
ties in the particle-hole channel are discouraged in normal 
circumstances with attractive microscopic interactions. 
Even if the interband couplings end up having repulsive 
character, the generic lattice and microscopic interaction 
potentials produce relatively small Ucv in comparison to 
Uc and Uv, so that a typical system with attractive in- 



teractions in d -t- e dimensions flows either to a charge- 
dynamics influenced insulator state, or toward particle- 
particle instabilities. With repulsive interactions, how- 
ever, the same kind of flows near the flxed points fea- 
turing Ucv take the system either to spin-dynamics influ- 
enced insulators, or toward the particle-hole instabilities. 

Finding the full structure of fixed points for any a > 
1/3 1 requires allowing all mixing interband couplings to be 
finite. Preliminary numerical calculations indeed reveal 
the existence of additional fixed points with finite mixing 
interactions and a < 1. However, a systematic search 
for these fixed points is very difficult due to the large 
parameter space and the highly non-linear nature of the 
RG equations that allow a to flow. The details of these 
flxed points are not crucial for the present discussion and 
will not be pursued further. 

Now we return to the possibility of band inversion 
which occurs for a < |/3|. First, we note that in nor- 
mal microscopic circumstances a is close to unity because 
the intraband couplings are larger than the interband 
ones C/™, n m. The RG flow further accentuates this 
situation as the flow of all is exponential. However, 
if the interband couplings are large enough in compar- 
ison to the intraband ones, we must cure the resulting 
effective band inversion by identifying the true low en- 
ergy quasiparticles, which must live at some new crystal 
wavevectors. The appropriate RG equations need to deal 
with more than two fcrmion flavors. Attractive interac- 
tions in such strong interband channels would naturally 
lead to paired states which spontaneously break transla- 
tional symmetry, while repulsive interactions would give 
rise to patterned exciton condensates. 

C. Transitions involving multiple fermion species, 
and universality classes 

The lowest energy quasiparticles in band insulators 
are often concentrated around multiple symmetry-related 
wavevectors in the Brillouin zone. For example, the sim- 
ple cubic periodic potential in three dimensions 

nr) = 2F[cos(^)+cos(|^)+cos(^) 

produces a band insulator with two fermions per site (for 
not too small V) whose lowest hole excitations live at 
fcv = (tt, 7r,7r)/aL in the valence band and lowest par- 
ticle excitations live at feci = (tt, 0,0)/aL and two other 
symmetry- related wavevectors fcc2, fec3 in the conduction 
band. An effective fermionic theory of this band insulator 
requires either one hole or three particle fields for generic 
pairing transitions of the type discussed in section III Al 
The discussion in section III Bl has to be extended to one 
hole and three particle fields in this case. 

An effective theory will generally include couplings 
among all of its fermion fields, and some of the couplings 
will have the same value by symmetries. As a proto- 
type theory we can take the action ^ allowing the la- 



10 



bels n,m . . . to identify any relevant fermion flavor. Like 
before, the RG analysis would reveal fixed points and 
run-away flows corresponding to same-flavor pairing and 
flavor-mixing instabilities. The latter kind could lead to 
supersolid phases in the particle-particle channel, or exci- 
ton condensates in the particle-hole channel, both bring- 
ing translational symmetry breaking and new universal- 
ity classes. On the other hand, the same-flavor pairing 
instabilities are the most likely outcome of generic attrac- 
tive interactions in lattice potentials due to the typically 
dominant couplings in the same-flavor channel. 

All superfluid transitions in the unitarity limit which 
involve only the same-flavor pairing always belong to the 
same universality class. This universality class can be 
characterized by critical ratios between pressure, temper- 
ature, energy per particle, and chemical potential (rela- 
tive to the band edge) at small but finite quasiparticle 
densities. A useful way of calculating the critical ra- 
tios involves applying a Hubbard-Stratonovich transfor- 
mation to the model ([Ij to decouple the short-range in- 
teraction, and then promoting the obtained two-channel 
model to an Sp(2iV) symmetry group by introducing N 
copies of the spinful fermion fields which couple to the 
same Hubbard-Stratonovich field. Fluctuation correc- 
tions to the mean-field thermodynamic functions take the 
form of 1/A^ expansions, so at least in the limit of large 
N one can obtain systematic perturbative expressions in 
the absence of a natural small parameter near unitarity. 
Taking the physical value = 1 and including only the 
lowest order correction ("Gaussian fluctuations" of the 
order parameter) already produces very good estimates 
in the uniform system^ldS. 

Provided that the inter-flavor scattering vanishes at 
the fixed point, trivial adjustments are needed to accom- 
modate multiple fermion flavors in the presence of a lat- 
tice, most notably in the ratios derived from extensive 
quantities, such as those containing pressure and energy 
density. For example, the critical pressure P at the finite- 
temperature T = Tc superfluid transition (in d = 3) 



/i < are universal functions of 



{P-Po)/N 



(2to)3/2T5/2„^ 



0.13188 



T=Ta 



0.4046 

A^ 



Oil/N"^) 



acquires a factor of Uf in the denominator on the 
left-hand side, the total number of low-energy particle 
and hole flavors in the Brillouin zone (Pq is the zero- 
temperature degeneracy pressure of the band insulator). 
The value of n/ depends on the bandgap Eg ; in the limit 
Tc ^ Eg both particle and hole flavors should be counted 
in Uf, otherwise only particles or holes are important 
based on the chemical potential. 

Another small adjustment of the uniform system 1/A^ 
expansions in Ref.'^'''*'' is needed in the critical ratios in- 
volving the chemical potential. We need to express the 
chemical potential /i relative to the nearest band edge. If 
the conduction band is nearest, then the finite quasipar- 
ticle density at zero temperature is obtained when > 0, 
while finite density requires /i < if the valence band is 
nearest. The critical temperatures at both fi > and 



T 



= 1.50448- 



2.785 

A^ 



This expression applies even in the limit Tc 3> Eg when 
both particles and holes are important, because this limit 
can be interpreted as ^ i?g. If the lattice depth is so 
small that the bandgap closes {Eg = 0), we must take the 
larger of the two values of obtained by measuring the 
chemical potential with respect to the overlapping "con- 
duction" and "valence" band edges. Additional phase 
transitions below Tc are possible for Tc ~ |^| » Eg, 
involving the onset of pairing in different channels: par- 
ticle, hole and interband, each characterized by its own 
order parameter (see previous section). Re-entrant be- 
havior can be anticipated in this regime when only one 
fermion species is paired at T = and another one is 
separated from the chemical potential by a gap smaller 
than Tc- Then, the thermal population of the fermions 
across the gap can lead to pairing in additional channels 
at < r^' < T < Tc. 

Anisotropy associated with low-energy quasiparticles 
at symmetry-transforming wavevectors in the Brillouin 
zone is equally easily treated. For any quasiparticle flavor 
with dispersion 



(21) 



we redefine momentum so that k'f/2m — kf /2mi, where 
m is a mass to be determined. The measure in path inte- 
grals acquires a factor of y/ (Hi rai)lm'^ from this change 
of variables, which can be absorbed into the redefinition 
of matter fields. This also leads to a renormalization of 
all interaction couplings. The choice 



(22) 



converts the quasiparticle dispersion into an isotropic one 
without renormalizing any fields or couplings. It is there- 
fore this geometric mean which should replace the mass 
in all 1/A^ expansions. 



III. DISCUSSION AND CONCLUSIONS 

We considered a band insulator subjected to pairing 
in the unitarity regime as a model system. The sim- 
plest realization of such a system is found in trapped 
neutral ultra-cold gases of alkali atoms placed in an opti- 
cal lattice. The density of atoms can be chosen to corre- 
spond to two atoms per lattice site in the central portion 
of the trap, while the strength of attractive interactions 
among them is routinely controlled by the Feshbach res- 
onance. A superfiuid transition from a thermally excited 
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band insulator has been already experimentally studied 
in this kind of a system in the vicinity of the BCS-BEC 
crossover—. 

The focus of our analysis was the characterization of 
the universal phase diagram featuring T = transitions 
between band insulators and superfluid states. In c? > 2 
dimensions we identified a BCS limit in which this tran- 
sition is pair-breaking, meaning that its universal prop- 
erties are transparently captured by a BCS-like theory. 
A special limiting case of the pair-breaking transition is 
found at unitarity, where all interaction effects become 
independent of microscopic scales, leading to the univer- 
sal dependence of critical temperature and other thermo- 
dynamic functions on the particle density in the super- 
fluid state. 

The BEC limit, found at any interaction strength in 
d = 2 or at sufhciently strong interactions in d > 2, 
brings a different universality class to superfluid tran- 
sitions. Fermionic excitations belong to high energies, 
so the effective theory capturing the transition has only 
bosonic fields. The transition occurs between the su- 
perfluid and a bosonic Mott insulator. The universality 
class is characterized either by the dynamical exponent 
z = 2 (generic bosonic mean-field transitions driven by 
the chemical potential), or z = 1 (XY transitions driven 
at fixed density of two fermions per lattice site) . 

The BCS and BEC limits considered here are rel- 
ative to a particular band insulator with a fixed lat- 
tice potential and particle density at zero temperature. 
While the particle density is finite in the ground-state, 
the unitarity regime between these BCS and BEC lim- 
its is found at zero quasiparticle density in the effective 
low-energy theory describing the band insulator. The 
full microscopic model includes short-range interactions 
and multiple fermion bands, with the chemical poten- 
tial residing in a bandgap. Integrating out high-energy 
fermions leaves behind the effective theory featuring at 
most two bands immediately adjacent to the chemical po- 
tential (the conduction and valence bands) . The remain- 
ing low energy fermions experience renormalized interac- 
tions, and may exist in multiple flavors as quasiparticles 
and quasiholes concentrated around different symmetry- 
related wavevectors in the first Brillouin zone (individu- 
ally having anisotropic dynamics). All of this complex- 
ity reduces to a few relevant interaction couplings in the 
vicinity of renormalization group (RG) fixed points that 
signify universal behavior, the most naturally occurring 
ones corresponding to unitarity in the same universality 
class as if the system were microscopically uniform. 

The physical meaning of these fixed points, revealed 
by RG, is the resonant scattering of quasiparticles. Mul- 
tiple flavors of quasiparticles give rise to multiple pos- 
sibilities for resonant scattering. An interesting discov- 
ered possibility is the resonant scattering between parti- 
cles and holes in the presence of repulsive interactions, 
the unitarity limit in the particle-hole channel separat- 
ing the regimes with non-existing and existing exciton 
bound states (excitonic "BCS" and "BEC" regimes re- 



spectively). Other possibilities not elaborated here also 
exist in generic circumstances with multiple fermion fla- 
vors, leading to translational symmetry breaking in or- 
dered states. However, most of these universal regimes 
may be inaccessible in realistic systems because they in- 
volve tuning either the details of lattice potentials, or 
short-range spatial features (at the lattice spacing scales) 
of the microscopic interaction potential. 

The notable exception are unitarity regimes in the uni- 
form particle-particle and hole-hole channels, which can 
be reached in cold atom systems using Feshbach reso- 
nances. The simplest to obtain is the unitarity at a tran- 
sition driven by the chemical potential, which is naturally 
found in a trapped gas of cold atoms at an interface be- 
tween the superfiuid and insulating atom clouds. In this 
case the RG identifies only one relevant interaction pa- 
rameter, which is tuned by the Feshbach resonance. The 
transitions driven at fixed density by changing the inter- 
action strength or lattice depth are harder to push to the 
full unitarity because there are two RG relevant oper- 
ators (particle-particle and hole-hole scattering lengths) 
which need to be tuned to their fixed point values. Nev- 
ertheless, manifestations of this kind of unitarity can be 
observed at finite temperatures if the critical tempera- 
ture is larger than the bandgap. 

The RG also provides an indication of the macroscopic 
properties of states away from the fixed points. If the 
strength of attractive interactions U is smaller by mag- 
nitude than its fixed-point value 1 17* | cx e in d -I- e dimen- 
sions, then a gapped fermion system is macroscopically 
a band insulator. Otherwise, the coupling U flows to- 
ward — oo under RG at finite length scales, implying the 
formation of Cooper pairs at short length scales before 
the onset of superfluidity at large length scales. It is in 
this manner that the fermionic RG predicts the existence 
of bosonic Mott insulators, but a bosonic effective the- 
ory is then required to access the superfluid transition at 
macroscopic scales. 

Perhaps the main significance of the presented RG 
analysis is that the most generic weak-coupling fixed 
points in fermionic theories, which control the univer- 
sal properties of insulating and superfluid phases, can be 
interpreted as resonant scattering. There is a distinction 
between appropriate "BCS" and "BEC" regimes in dif- 
ferent kinds of pairing channels, in terms of the existence 
of appropriate two-quasiparticle bound states. The run- 
away flows of interaction couplings in the "BEC" regimes 
signify the emergence of correlated insulating states sep- 
arated from ordered phases by transitions in bosonic uni- 
versality classes. In some circumstances these "BEC- 
limit" insulators may be thermodynamic phases, such as 
a valence-bond crystal or a spin liquid adjacent to an an- 
tiferromagnet (condensate of excitonic "molecules"), or 
a charge-density wave adjacent to a superconductor. 

Therefore, the presented model and analysis provide 
a direct insight into the possibilities for the develop- 
ment of strong pairing correlations in fermionic lattice 
systems. The emergence of boson-dominated superfluid 
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transitions among fermions and the corresponding uni- 
versality classes can be traced down to the well known 
physics of BCS-BEC crossovers. Even if interactions are 
not strong enough to bring the system close to its uni- 
tarity limit in empty space, the presence of a lattice frus- 
trates the motion of particles and promotes interaction 
effects, effectively pushing the system toward its lattice 
unitarity^^'^^. Furthermore, in two dimensions there is 
no BCS limit strictly speaking. Two quasiparticles in- 
jected into the conduction band will form a bound state 
regardless of how weak the attractive interactions are. Of 
course, the size of this "vacuum" bound state might be 
much larger than the spacing between particles, but this 
does not preclude the bosonic universality of the super- 
fluid transition. 

One potentially important aspect of this is that a 
conceptually similar situation is found in cuprate high 
temperature superconductors. Cuprates are quasi two- 
dimensional systems in which the underdoped normal 
state (pseudogap) exhibits gapped fcrmionic quasipar- 



ticles, albeit with a specific rf-wave pairing symmetry 
and a gap of completely different origin than in this pa- 
per. A number of unconventional properties of cuprates 
can be qualitatively understood as being related to a 
fluctuation-driven transition. 
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